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ABSTRACT 

We combine lensing, stellar kinematic and mass-to-light ratio constraints to build a two- 
component (luminous plus dark) mass model of the early-type lens galaxy in PG1 115+080. 



We find a total mass density profile steeper than r , effectively p ex 



with y 



2.35 ± 0.1 ± 0.05 (random + systematic). The stellar mass fraction is /» = 0.67Z U Z 2 U 5 ± °- 03 
inside the Einstein Radius (Re«1.2 effective radii). The dynamical mass model breaks the de- 
generacies in the mass profile of the lens galaxy and allows us to obtain a value of the Hubble 
Constant that is no longer dominated by systematic errors: H = 59±7 2 ± 3 kms -1 Mpc -1 
(68% C.L.; il m = 0.3, O A =0.7). The offset of PG1 115+080 from the Fundamental Plane 
might indicate deviations from homology of the mass profile of some early-type galaxies. 
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1 INTRODUCTION 

Time delays between multiple images of gravitational lenses pro- 
vide the opportunity to measure the Hubble Constant Ho (Refsdal 
1964) independently of local distance scale methods which rely on 
uncertain calibrations (Saha et al. 2001; Freedman et al. 2001). Al- 
though the method is attractive for its reliance on general relativity 
alone and being a one-step global measurement, degeneracies in- 
herent to the mass distribution of the lenses seem hard to break 
without external information (e.g. Koopmans 2001). 

We recently showed that a combination of stellar kinematics 
and gravitational lensing can be used to place tight constraints on 
the mass distribution of early-type (E/S0) lens galaxies inside their 
Einstein radius (Koopmans & Treu 2002, Treu & Koopmans 2002; 
KT02 and TK02). In particular, it was found that the total mass den- 
sity profiles of the lens galaxies in MG2016+1 12 and 0047-281 are 
very close to p oc r -2 . Both lens systems lie on the Fundamental 
Plane (FP; Dressier et al. 1987; Djorgovski & Davis 1987) at their 
respective redshifts, and are part of the sample that we selected for 
the Lenses Structure and Dynamics (LSD) Survey, currently being 
carried out at the Keck Telescope. 

PG1115+080 (,z=0.31; Weynman et al. 1980) is not part of the 
LSD Survey - which focuses mainly on relatively isolated early- 
type galaxies - because of the presence of a massive compact group 
nearby that could affect the mass distribution of the lens galaxy 
through interaction. However, this system is particularly interesting 
because time delays between the lensed quasar images (z s — 1.71) 
are available (Schechter et al. 1997; Barkana 1997), which can be 
used to measure Ho. The geometry of this lens system has been 
modelled in detail (Schechter et al. 1997; Keeton & Kochanek 



1997; Courbin et al. 1997; Saha & Williams 1997; Impey et al. 
1998, hereafter 198; Williams & Saha 2000; Kochanek, Keeton & 
McLeod 2001; Saha & Williams 2001; Zaho & Pronk 2001), illus- 
trating how strongly the value of Ho depends on the mass profile 
of the lens galaxy. Whereas par - mass models yield values of 
Ho < 50 km s -1 Mpc -1 , steeper or constant M/L mass models can 
yield values up to 60-70 km s -1 Mpc -1 . 

In this Letter, we apply the methods developed for the LSD 
Survey to the lens E/S0 galaxy in PG1 1 15+080 to break the mass- 
profile degeneracy and perform an accurate measurement of the 
Hubble Constant, minimising this dominant source of systematic 
uncertainty. Based on the lensing and kinematic constraints (i.e. 
the published velocity dispersion measurement by Tonry 1998), we 
build a mass model (§2) which we then use to determine Ho (§3). 
We summarise and discuss our results in §4. Throughout, we de- 
fine h = Ho/100 kms -1 Mpc -1 , and we assume f2 m = 0.3 and 
JIa = 0.7. The cosmography affects Ho only at a few percent level. 



2 THE MASS DISTRIBUTION 

2.1 Spectrophotometric properties of the lens 

Hubble Space Telescope (HST) images of the gravitational lens 
system PG1 115+080 are available from the HST archive. We se- 
lect the images with better resolution and sampling in each of the 
available bands: 4 x 640s exposure with the Near Infrared Camera 
and Multi Object Spectrograph (NICMOS; GO-7496; PI: Falco) 
Camera 2 (NIC2) through filter F160W (198); 10 exposures on the 
Planetary Camera (PC) of the Wide Field and Planetary Camera 2 
(WFPC2) through filter F814W for a total exposure time of 15060 s 
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Table 1. Surface photometry of the lens galaxy. Rest frame quantities 
through filters B and V are computed for h = 1, f2 m = 0.3, Qjy = 0.7 
and E(B-V)=0.041 (Schlegel et al. 1998) 





F160W 


F814W 


m (mag) 


16.36 ±0.15 


18.47 ±0.10 


SB e (mag arcsec -2 ) 


17.75 ±0.25 


20.11 ±0.10 


R e (arcsec) 


Of/76 ±0'.' 12 


0'.'85 ± 0"07 




V 


B 


M (mag) 


-20.89 ±0.12 


-20.05 ±0.12 


SB e mag arcsec ~~ 2 


19.83 ±0.12 


20.67 ±0.12 



(HST GO-6555; PI: Schechter); 8 x 400s exposures on the WFPC2- 
PC through filter F555W (HST GO-7495; PI: Falco) 

The images were first reduced using a series of IRAF tasks 
based on the DRIZZLE package (Fruchter & Hook 2002). This pro- 
cess yielded maps of bad pixels, cosmic rays, saturated pixels, and 
pixels affected by horizontal smearing for each individual expo- 
sure. Then, we used a series of IDL scripts to subtract the QSO 
multiple images from each individual exposure, after masking the 
deviant pixels. For each exposure we used a library of subsampled 
synthetic Point Spread Functions (PSFs) computed with TINY TIM 
6.0 (Krist & Hook 2001) covering a range of the relevant parame- 
ters, such as jitter and focus offset to simulate HST breathing. The 
PSF giving the smallest residuals was used to remove the QSO im- 
ages from each exposure. The individual exposures were finally 
combined to produce an image of the lens galaxy. The lens galaxy 
is imaged at high signal to noise ratio in the F160W and F814W 
images, which we used to measure its structural parameters. The 
structural parameters are listed in Table 1. An almost complete 
ring is visible in the F160W image (see 198). Unfortunately, the 
signal-to-noise ratio is significantly lower in the F555W exposure 
and structural parameters could not be reliably measured. However, 
the colour F555W-F814W=1.82 ± 0.02 can be reliably measured 
within a fixed aperture of radius 0'.'6. 

Rest frame quantities were obtained as described in Treu 
et al. (2001b) using the F814W photometry, the F555W-F814W 
colour, and correcting for galactic extinction with E(B-V)=0.041 
(Schlegel, Finkbeiner & Davis 1998). We use the parameters in the 
F814W filter because it is much closer in wavelength than F160W 
to the standard rest frame B and V bands, and the effects of possible 
colour gradients are therefore minimised. 

Using the stellar velocity dispersion (Tonry 1998) of 
281 ±25 kms -1 measured inside a l'.'O squared aperture, which 
corresponds to a central velocity dispersion a — 293 ± 26 km s _1 
in a circular aperture of radius r e /8, we place the lens galaxy of 
PG1 115+080 in the FP space (i.e. the space with axis loger, SB e , 
and log R e ) and compare the position of the galaxy to the FP as 
defined by the largest sample of galaxies at comparable redshift, 30 
E+S0 galaxies in the cluster CL 135 8+62 at z = 0.33 (Kelson et 
al. 2000). As shown in the left panel of Fig. 1, the lens galaxy is 
dimmer than cluster galaxies with the same velocity dispersion and 
effective radius (see also Kochanek et al. 2000). The lens galaxy 
is dimmer even than the galaxies in the Coma Cluster (Bender et 
al. 1998) at variance with the observed trend that the luminosity of 
E/S0 galaxies at a given radius and velocity dispersion increases 
as a function of increasing redshift (e.g. Treu et al. 2002). In the 
right panel of Fig 1 we show the residuals from the FP of Coma 
for the lens galaxy and the galaxies in cluster CL1358+62. The 
lens galaxy lies on the dimmest tail of the distribution, and is offset 
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Figure 1. Left: Comparison between the properties of the lensed galaxy in 
PG1115+080 (filled square) and the FP of CL1358+62 (empty pentagons). 
The hatched square is plotted using ctsie = 219 ± 5kms _1 as central 
velocity dispersion. Right: Offset from the FP of Coma (zero), in the rest 
frame B band for the galaxies in CL1358+62 and PG1 115. Open histogram 
represents galaxies in CL1358+62, solid and hatched histogram PG1 1 15 as 
in the left panel. 



by 4-5 times the rms from the median/average CL1358+62 value. 
Since nothing in the morphology or colours (e.g. the lens has B- 
V=0.84 ± 0.05, whereas the average colour of the galaxies in clus- 
ter CL1358 is B-V=0.86) of the lens galaxy seems to indicate the 
presence of dust or anomalous stellar populations, we will assume 
that the offset arises from structural differences. 

As in TK02 and KT02 we use the evolution of the effec- 
tive mass-to-light ratio A log(M/L_g), inferred from the local FP 
studies, to estimate the stellar mass to light ratio M*/Lb of the 
lens galaxy. Assuming that \og(M r /L B )z = log(M*/L s )o + 
Alog(M/L_g) (see discussion in Treu et al. 2001a; Treu et al. 
2002; KT02), we use the average local value 12.0±4.1 ft M /L s , 
(Gerhard et al. 2001 and references therein) and the average evolu- 
tion observed for field E/S0 galaxies A log M/Lb = -0.22 ± 
0.04 at 2=0.31 (Treu et al. 2002) to find M,/L B = 7.2 ± 
2.5ft M /Ls, for the lens galaxy. The observed evolution 
for cluster galaxies by van Dokkum et al. (1998) yields 8.5 ± 
2.9 ft Mq/Lb_0. In the context of our hypothesis of a structural 
origin of the FP offset, we use the average observed values and 
not the evolution derived from the galaxy itself to estimate its stel- 
lar mass to light ratio (which would yield M,/Lb = 14.6 ± 
5.1ftM /L B , ). 



2.2 Dynamical and lensing models 

We calculate the expected stellar velocity dispersion of the lens 
galaxy in PG1 1 15+080 - to compare then with the observed value, 
§2.1) - by modelling its mass distribution with two spherical com- 
ponents 1 , one for the luminous stellar matter and one for the dark- 
matter halo. The luminous mass density profile is described either 
by a Hernquist (1990; HQ) model or a by Jaffe (1983; JF) model. 
The dark-matter distribution is modeled with a generalised ver- 
sion of the Navarro, Frenk, & White (1997; NFW) density profile. 
CDM simulations show that the break radius (r&) is much larger 



1 See, e.g., Saglia, Berlin & Stiavelli (1992 and Kronawitter et al. 200) for 
discussions on the accuracy of the spherical approximation. 
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Table 2. NICMOS astrometry and relative fluxes. Plate scales are 
O'.'0760261/pix and O'.'0753431/pix along x and y respectively, taken from 
the NICMOS history tool available at STScI web site. The position angle of 
the NICMOS image is 68.7594 degrees. 



Object 


ARA 


A DEC 


S h/ S h 


Al 


+1"328 ± 0'.'006 


-2" 042 ± 0"006 


2.29±().()8 


A2 


+1//472 ± 0'.'005 


-1"581 ± r/.'005 


3.44±0.15 


B 


-tf'338 ± 0'.'006 


-1"965 ± 0'.'005 


0.76±0.04 


C 


+tf'000 ± 0'.'006 


+(y.'000 ± 0'.'004 


1.00 


G 


-0'.'381 ± 0'.'007 


-1"345 ± r//007 





than the Einstein Radius (e.g. Bullock et al. 2001). Hence, effec- 
tively pd(r) oc r~ 7 inside the Einstein radius, where 7 is the 
inner slope of the dark matter halo (7=1 for NFW). We further 
adopt an Osipkov-Merritt (OM; Osipkov 1979; Merritt 1985a,b) 
parametrisation of the anisotropy of the stellar velocity ellipsoid, 
P(r) = 1 — a 2 /a 2 = r 2 /(r 2 + r 2 ), where ae and a r are the 
tangential and radial components of the stellar velocity dispersion 
and n indicates the anisotropy radius (see KT02 for discussion). 
The line-of-sight velocity dispersion is determined by solving the 
spherical Jeans equation, correcting for the average seeing of (X'8, 
and averaging the velocity dispersion profile - weighted by the sur- 
face brightness - inside the 1 □" spectroscopic aperture. 

A very similar mass model is used to reproduce the lensing 
geometry with two minor differences: (i) the lens model allows for 
ellipticity of the dark-matter halo (b/a > 0.9 in all cases) whereas 
the luminous mass profile is spherical (consistent with the obser- 
vations; §2.1; 198); (ii) the luminous mass profile is modelled with 
a HQ or pseudo-Jaffe profile (Keeton 2001). The latter is analyti- 
cally tractable and differs only marginally from the JF profile. An 
additional component, accounting for the nearby massive compact 
group, is necessary to reproduce accurately the lensing geometry 
(e.g. Schechter et al 1997). We model the group as an isothermal 
sphere (i.e. p oc r~ 2 ) or a NFW mass distribution. The break ra- 
dius of the NFW profile is chosen to be 10", consistent with results 
from CDM simulations (e.g. Bullock et al. 2001) for the observed 
velocity dispersion of the group (~330 km s _1 ; Tonry 1998). 

2.3 Constraining the mass profile of the lens 

We now use the available observations to constrain the free param- 
eters of the models introduced in §2.2, with the overarching goal 
of determining accurately the mass profile and hence the value of 
Ho. To this aim we use four sets of observables. (i) The relative as- 
trometry of the multiple images and lens galaxy, derived from the 
NICMOS images using PSF fitting, isophotal fitting, and centroid- 
ing (Table 2; names are as in 198). The uncertainties include the sta- 
tistical component and the systematic difference between the three 
techniques. Our astrometry is consistent within the errors with that 
obtained by 198. (ii) The observations described in §2.1, namely 
the effective radius, total luminosity, velocity dispersion and stel- 
lar mass-to-light ratio of the lens galaxy, (iii) The time delay be- 
tween images B and C Atsc = 25.0 ±1.6 days (Barkana 1997; 
Schechter et al. 1997). (iv) The relative fluxes of the images with 
20% errors (Table 2; see 198). The precise choice of the errors on 
the flux-ratios are less important, because the more precise image 
positions are the dominant contraints on the lens mass model. 

Specifically, we use these observations and the dynamical 
model to constrain the slope of dark matter halo (7) and the frac- 
tion of stellar mass within the Einstein Radius (/*). The relative 



astrometry yields the projected mass of the lens galaxy within the 
Einstein Radius - corrected for the projected mass contribution of 
the group - which is virtually independent of the mass profile of the 
galaxy or that of the nearby group for the two group mass profiles 
that we adopt (see §2.2). 

We find a mass of M(R < R E = l'/04) = (1.19 ± 0.06) x 
10 11 /i -1 M0 - corresponding to ctsie — 219 ± 5 kms -1 for 
a singular isothermal ellipsoid - consistent with results from 198. 
The Einstein radius Re=1'.'04 corresponds to Re = 3.32 ft -1 kpc 
at the redshift of the lens. The half-light radius of the luminous 
component is set equal to the effective radius of the best fitting 
R 1/A model (§2.1). For any given value of the anisotropy radius n, 
we determine the likelihood of each (/„ ,7) pair, by comparing the 
observed stellar velocity dispersion with that from the dynamical 
model, and assuming Gaussian error distributions. In the limit of 
negligible stellar mass fraction (/„ — > 0), the "effective slope" (see 
TK02 or KT02) of the lens-galaxy mass profile is 7' = 2.38 ± 0.08 
and 7 ' = 2.31±0.08 (68% CL) for the isotropic (n = oo)HQand 
JF luminous profiles, respectively. These values are lower by 1-2% 
for non-isotropic models with n — R c . This suggests that a steep 
total mass profile is required for the lens galaxy to explain the large 
observed stellar velocity dispersion and is expected to increase Ho 
by 30^-0% (Wucknitz 2002). For completeness, we note that a total 
mass profile with 7' = 2 can be excluded at the 99.9%(99%) CL 
for both the isotropic(non-isotropic) JF and HQ models. 

If we now include the constraints on the stellar (cluster) mass- 
to-light ratio (§2.1; see TK02 or KT02 for details) - which gives 
the stellar mass fraction inside the Einstein radius independent of 
Ho - we find dark-matter slopes of 7=2. 35 ± 0.25 (68% CL) and 
7=2.60 ± 0.20 for the isotropic JF (shown as an example in Fig. 2) 
and HQ profiles, respectively. Mildly anisotropic JF models with 
ri = R e lead to 7 = 2.20! ( l ;g \ whereas similar HQ models 
give 7 = 2.50l^. The stellar mass fraction for a JF profile is 
/„ = 0.70lo;2° (68% CL; compare, e.g., the results in Bertola et 
al. 1993 for local E/S0) inside the Einstein radius (0.05 lower for 
HQ profile). This fraction lowers by 0. 1 when we use the field M/L 
evolution (§2.1). 

For the HQ profile with its relatively shallow inner luminos- 
ity density profile, the dark-matter halo density dominates at lens- 
galaxy radii < 0.5 kpc (in fact, this happens even for HQ profiles 
in an r -2 total mass profile). For this reason we regard this a less 
likely mass model, although we will continue to use it to illustrate 
the effect of a range of inner luminosity-density profiles on the de- 
termination of H . A small core radius of ~ 0.01 kpc for the steeper 
JF models is sufficient to avoid this problem and has only negligi- 
ble effect on the stellar dispersion. We also find that r~ 2 total mass 
profiles with neither a constant /3(r) (negative or positive, i.e. tan- 
gentially or radially anisotropic) nor an Osipkov-Merritt model at 
the limit of radial instability (e.g. Nipoti, Londrillo & Ciotti 2002) 
can reproduce the observed stellar velocity dispersion. 

3 THE HUBBLE CONSTANT 

Having the constraints on /* and 7 in hand - we can tightly con- 
strain the lens models and determine Ho (§3), thereby including in 
its error budget the uncertainties on the mass profile of the lens. 
We use the lens-code from Keeton (2001) to model the system with 
the data described above. As a test, we are able to recover the re- 
sults from 198 within the errors with both their astrometry and ours. 
The stellar mass fraction and the slope of the dark-matter halo in- 
side the Einstein radius are set equal to those determined in the 
previous paragraph. This leaves as free parameters: the ellipticity, 
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Table 3. Values of H from PG1 115+080. For each combination of 
mass models for the lens galaxy and the nearby group we list the value 
of Ho in kms~ 1 Mpc~ 1 along with its random uncertainty (JF=Jaffe, 
HQ=Hernquist, ISO=isotropic, OM=Osipkov-Merrit anisotropy with r-j = 
R e , SIS=singular isothermal sphere, NFW=Navarro, Frenk & White). All 
models have \ 2 in the range 2.8-2.9 (without the ring constraints). 







JF-ISO 


HQ-ISO 


JF-OM 


HQ-OM 


group 
group 


SIS 
NFW 


57+J, 2 
59t* 2 


61+1 2 

621" 


56+11 
57ll 2 


60+12 
62+12 



the position angle, and mass inside the Einstein Radius of the dark- 
matter halo; the source position; the position and mass of the nearby 
group. For each mass model, we find the set of free parameters that 
minimises the x 2 - From the resulting minimum-x 2 lens model, we 
determine Ho by comparing the model time-delay between images 
B and C with the observed time delay (§3). We stress that all con- 
straints on the lens model are independent of Ho. 

Additional constraints can be obtained by using the shape of 
the Einstein ring (Kochanek et al. 2001), as determined by trac- 
ing its peak surface brightness in 47 independent points. We do 
not use the brightness maxima or minima of the ring, which are 
difficult to measure accurately. We assume an error of 0'.'02 on 
the ring peak-brightness positions, chosen to yield a minimum- 
X 2 approximately equal to the number of degrees of freedom for 
the r~ 2 lens model. In this way, we find minimum-x 2 values of 
56, 55 and 55 for models where the lens galaxy is modeled as 
a single singular isothermal ellipsoid (198), a single Pseudo-JF 
(i.e. constant M/L) profile with half-light radius equal to R e (i.e. 
a = (4/3) R e = 1'.'13; see Keeton 2001), and our best HQ model 
embedded in a steep dark-matter halo (see §2.2), respectively (a sin- 
gular isothermal group is assumed). Remarkably, the x values are 
indistinguishable and the values of Ho in all three cases are found 
to be similar to those without the ring as an additional constraint. 
We therefore conclude that the Einstein ring does not allow us to 
distinguish between an r~ 2 mass model, a constant M/L pseudo- 
JF model with break radius a < 2'.'0, or our best dynamical mass 
model, nor does it improve significantly the accuracy on Ho (see 
also Saha & Williams 2001). Finally, we note that all our models 
giverABC = AtAc/AteA ~ 1.3, consistent with Barkana (1997). 

Tab. 3 lists the values of Ho obtained for a variety of mass 
models (HQ to JF luminous mass profiles, isothermal to NFW 
group mass profiles, isotropic to non-isotropic stellar mass distri- 
butions), all consistent with the available observations (see §2.3). 
This relatively large set of models leads to a tight range of Ho=56- 
62kms _1 Mpc _1 . 

To compute the random errors we include the contributions 
from uncertainty in the images and lens positions, in the time-delay, 
in the image flux-ratios, in the slope of the dark-matter halo and in 
the fraction of stellar mass inside the Einstein radius. The errors 
from the lens model and dynamical model are determined sepa- 
rately, assuming they are independent. In all cases the errors on Ho 
due to the lens models significantly dominate those due to the dy- 
namical models. In particular, the uncertainties on the galaxy po- 
sition and time-delay dominate the error budget on Ho. We com- 
bine the random errors from the lensing and dynamical models and 
list the 68% CL limits in Tab. 3. Because we re-optimize x 2 until 
Ax 2 = 1, the error-ranges include the effects of any parameter de- 
generacies (i.e. internal versus external shear) within the context of 
our mass model. 

We adopt, as best estimate of H from PG1115, the average 
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Figure 2. Likelihood contours of fraction of stellar mass within the Ein- 
stein radius vs. dark matter slope (7) for the isotropic (left panels) and 
non-isotropic (right panels) JF model. We note that both the slope of the 
dark-matter halo (7) and the stellar mass fraction enclosed by the Einstein 
radius are independent from Ho. 

value from the eight models considered in Table 3, for the ran- 
dom uncertainty the average of the random uncertainties and for 
the systematic uncertainty the semi-difference between the maxi- 
mum and minimum values, i.e. H = 59+ 12 ± 3 kms" 1 Mpc^ 
(random/systematic) . 

4 SUMMARY & DISCUSSION 

We have combined mass-to-light ratio constraints with lensing and 
dynamical analyses - using the stellar velocity dispersion measure- 
ment from Tonry (1998) - to constrain the radial mass profile of 
the lens galaxy of PG1 1 15+080 and break its degeneracy with the 
value of Ho, as determined from the measured time delay. Our main 
results are: 

(1) A steep total mass profile with pocr~ 7 , where 
7'=2.35±0.1±0.05, and an isotropic or mildly radial velocity el- 
lipsoid, successfully reproduces the stellar velocity dispersion and 
the lensing constraints. A two-component mass model, with a lu- 
minous component with stellar M*/Lb as determined from the 
evolution of the FP, satisfies all the constraints only for a relatively 
steep (compare TK02 and KT02) inner slope of the dark matter 
halo and a considerable fraction of stellar mass within the Einstein 
Radius, /„ = 0.67to' 2 5 =t 0.03. This model explains the offset 
of PG 11 15+080 from the FP in terms of structural differences and 
does not require strong radial anisotropy of the stellar velocity dis- 
tribution. 

In previous work, ad hoc distributions of stellar orbits have 
been constructed inside a logarithmic potential also leading to a 
high stellar velocity dispersion (Romanowsky & Kochanek 1999). 
However, if the enclosed mass is kept fixed (given by the lens 
model), these orbits need to be extremely radial and instability 
issues may arise. Recent numerical analyses show that if E/S0 
galaxies are offset from the FP through radial anisotropy by more 
than the observed scatter in the FP, they become unstable and 
evolve back to the FP (Nipoti et al. 2002). Tangential orbits always 
lead to a lower central velocity dispersion. Our mass model ex- 
plains the offset by breaking the homology in the mass distribution, 
and is consistent with the anisotropy structure expected in galaxy- 
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formation scenarios, with close to isotropic orbits in the central re- 
gions (e.g. van Albada 1982; see also the observations in Gerhard 
etal. 2001). 

Our similar analyses of two other lenses (MG2016+112 and 
0047-281; TK02 and KT02) indicate that both lie on the FP of 
field E/SO galaxies and not only have total mass profiles with slope 
—2 within 5%, but in the case of 0047-281, a spatially resolved ve- 
locity dispersion profile allows us to rule out tangential anisotropy 
or significant radial anisotropy in the center of the galaxy. Data for 
other lenses with quality comparable to 0047-281 are being col- 
lected by the LSD Survey, to provide more information on the mass 
distribution and orbital structure of E/S0 (lens) galaxies. 

(2) Applying these dynamical constraints to the lens mass 
models of PG1115+080, we find H = 591^ ± 3 kms" 1 Mpc" 1 
including random (68% CL) and systematic uncertainties. To 
achieve higher accuracy it will be necessary to further reduce the 
dominant sources of uncertainty, i.e. errors on astrometry and time 
delay. Note that this value is significantly higher than the value of 
Ho = 44±4 kms _1 Mpc _1 found by 198, who assumed an r~ 2 
mass profile. 

Comparing to other lensing measurements of Ho, our value 
is consistent with those derived from B 1608+656 (Fassnacht et 
al. 1999; Koopmans & Fassnacht 1999; Fassnacht et al. 2002) 
and B 1600+434 (Koopmans et al. 2000; Kochanek 2002), both 
assuming r~ 2 mass profiles. However, in other cases values of 
Ho < 50km s -1 Mpc -1 have been found using r~ 2 mass profiles 
as well (e.g. Kochanek 2002). We therefore conclude that the r~ 2 
approximation might not always appropriate for the determination 
of Ho from time-delays, given the intrinsic scatter in the total mass 
profile and fraction of dark matter in the central 1-2 R e of E/S0 
galaxies (see e.g., Bertin et al. 1994; Carollo et al. 1995; Gerhard 
et al. 2001, for ranges of mass profiles of E/S0 galaxies), particu- 
larly if there is evidence for offsets from the FP. We propose that 
the apparent differences between Ho, inferred from different lens 
galaxies modelled with the same mass profiles, is therefore the re- 
sult of structural, not kinematic, non-homology of E/S0 galaxies, 
which might be related to differences in their formation or interac- 
tion histories (e.g. field versus group/cluster). Hence it appears that 
additional constraints, such as stellar kinematics, are essential to 
precisely determine the mass profile - hence Ho - for at least some 
lens-systems with measured time-delays. 

Finally, as compared to local determinations of H , our value 
is consistent with Ho = 59 ± 6 kms -1 Mpc -1 from Saha et al. 
(2001; and references therein), but only marginally consistent with 
the final value, H = 72 ± 8 kms -1 Mpc~\ from the HST Key- 
project (Freedman et al. 2001). 
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